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ABSTRACT 

The non-uniform longitudinal flow model (NUFM) propsed recently is 
extended to include also transverse flow. The resulting longitudinally non- 
uniform collective expansion model (NUCEM) is applied to the calculation 
H ' of rapidity distribution of kaons, lambdas and protons in relativistic heavy 

ion collisions at CERN-SPS energies. The model results are compared with 
the 200 A GeV/c S-S and 158 A GeV/c Pb-Pb colhsion data. The central 
dips observed in experiments are reproduced in a natural way. 



PACS number(s): 25.75.-q 



^Work supported in part by the NSFC under project 19775018. 



I. Introduction 

Relativistic heavy-ion collision offers a unique opportunity to study the hot and dense 
matter under controlled laboratory conditions |I|. Hadronic spectra from these re- 
actions reflect the dynamic of the hot and dense zone formed in the collision. The 
baryon density, established early in the reaction, is an important factor governing the 
evolution of the system. Comparison of model predictions with measured rapidity and 
transverse momentum distributions constrains the possible dynamical scenarios of the 
reaction, such as those for longitudinal and transverse flow. In addition, the mechanism 
by which the incoming nucleons lose momentum during collision (baryon stopping) is 
an important theoretical problem 0. 

The rich physics of longitudinal and transverse flows is due to their sensitivity to the 
system evolution at early time. Due to the high pressure of the system created during 
a heavy-ion collision, particles might be boosted in the transverse and longitudinal 
directions. The expansion and cooling of the heated and highly compressed matter 
could lead to a considerable collectivity in the final state. The collective expansion 
implies space-momentum correlation in particle distributions at freeze-out. 

In this paper we present a study of the experimental results from the central colli- 
sions of Pb-Pb at 158 A GeV/c and S-S at 200 A GeV/c measured by NA49 and NA35 
at CERN SPS 0|0. The primary goal of the study is to provide a chance to compare 
the stopping power and radial symmetrical flow in collision systems of different sizes. 
Note that at SPS energies the duration time of the evolving system will be longer in 
comparison with that at the AGS energies, and it is even more reasonable to assume 
that a thermalized system expands collectively both in the longitudinal and at the 
same time in the transverse direction. 

A big challenge for thermal models coming from experimental data is the central 
dips in the rapidity distributions of heavier particles observed in the central collisions 
of nuclei both at AGS and at SPS energies. It is well known that models with uniform 



longitudinal flow [|12l cannot account for such dips. In order to reproduce the experi- 
mentally observed central dips, a parameterization of baryon chemical potential /ij as 
function of flow rapidity rj is introduced in Ref. [^]. However, the physical meaning of 
such a parametrization is unclear. 

In the present paper we will construct a longitudinally non-uniform collective ex- 
pansion model (NUCEM) basing on the physical argument that the flreballs produced 
in the nuclear collision will keep some memory on the motion of the incident nuclei. 
The non-uniform longitudinal flow will be parametrized using a geometrical picture 
proposed in Ref. |0|. The results of model calculation are compared with the exper- 
imental data on the rapidity distributions of baryons and strange particles in 200 A 



GeV/c S-S and 158 A GeV/c Pb-Pb collisions. 

The longitudinally non-uniform collective expansion model will be formulated in 
section II. The results of model calculation are presented and compared with the ex- 
perimental data in section III. A short summary and conclusion will be given in section 
IV. 

II. Longitudinally non-uniform collective expansion flow 

In the non-unifrom longitudinal flow model (NUFM) proposed in Ref only longitu- 
dinal flow was taken into account. In fact, the one dimensional expansion presumably 
dominates initially because of anisotropic initial conditions. But it is inconsistent to 
assume that a thermalized system expands collectively only in longitudinal direction 
without generating transverse flow from the high pressure in the hydrodynamic system. 
In case of central collisions the hydrodynamic flow should be computed in at least (2+1) 
dimensions if azimuthal symmetry is maintained. Beside these theoretical necessities 
for the inclusion of transverse flow [R|, there are also evidence for the existence of this 



flow from the experimental phenomena P|||10||. 



For the description of particle production in the collective expansion of flre-ball we 
start from the formalism of Cooper and Frye [|1^ which describes the single-particle 
spectrum as an integral over a freeze-out hypersurface, thus summing the contributions 
from all space-time points at which the particles decouple from the flreball: 

where g is the degeneracy factor and f{x,p) the momentum distribution at space-time 
point X. In thermal models one takes f{x,p) as a thermal equilibrium distribution and 
determines Uf by a freeze-out criterion for thermal decoupling [|1^. At freeze-out, the 



Boltzmann approximation is sufficient, but we allow for a space-time dependence of 
the temperature T, the chemical potential /i, and the flow velocity u'^: 

/(x,p)=exp(^ T{x) )■ ^'^ 

Focusing on central collisions, we assume azimuthal symmetry of the spatial geometry 
and momentum distributions. 

The assumption of longitudinal Bjorken flow [13] suggests to use longitudinal proper 
time r = yjt"^ — z^ and space-time rapidity rij\ = tanh~^(2;/t) as suitable variables in 
the t-z plane. The transverse radial coordinate is denoted by r. In an azimuthally 
symmetric geometry of this kind it is practical to decompose the velocity fleld in the 



following way: 

u'^ = (cosh r^t cosh r^i, erSinh?7t, cosh ?7t sinh?7r)- (3) 

Here ej. is the 2-dimensional unit vector in the radial direction, rji = rii{t,r,z) is the 
longitudinal flow rapidity by which each volume element on the z-axis moves relative 
to the center of mass, and rjt = Tit{t, r, z) is the rapidity corresponding to the transverse 
flow of the volume element at position (r, z) as seen from a reference point at z on the 
beam axis moving with the local flow velocity there. 

The momentum of the particle in the center-of-fireball system can be parametrized 
as follows: 

u^Pi, = rrit cosh(t/ — rji) cosh rjt — pt sinh rji cos(0 — ip) (4) 

Because of azimuthal symmetry in central collisions we can integrate over making 
use of the modified Bessel function Io{a) = (27r)^^ / e°"^°^'^d(f). 



The geometry of the freeze-out hypersurface Cf is fixed as follows: In the time 
direction we take a surface of constant proper time, r = rf. In t^i direction the freeze- 
out volume extends only to a maximum space-time rapidity rjo, which is required by 
the finite available total energy and breaks longitudinal boost-invariance propsed by 
Bjorken. In the transverse direction the boundary is given by -Rf, which describes a 
cylindrical fireball in the rj-r space. 

Having specified the freeze-out geometry and the distribution function at freeze-out 
we obtain the following thermal single-particle spectrum: 



d'^N g 



-mtTf 



- ,..,,, /dr^irdre'^/^e-"^"^^^^-^") cosh (y-r7i)^o(«), (5) 

where a = {rrit/T) cosh rjt, a = (pt/T) sinh rjt. After integrating over rrii, we get the 
rapidity distribution as follows: 

^ I dml I d'qi /'drrf(r)re^/^e-"=°^^(^-''') 



dy 



-vo u 

xmtCosh(y-?7i)Jo(«). (6) 



Here m}.°, m^^ denote the experimental limits in which the spectrum was measured. 
The freeze-out radius -Rf and the longitudinal extent of the fireball is fixed via the 
finite interval (—r]Q,r]Q). rjt = tanh~ Pt is the rapidity of transverse fiow. Replacing 
the transverse fiow velocity profile by its radial average, we get 

dN R2 ""'" ''° 

oiy^gniM j ^^2 ji. / dr7cosh(|/-77i)e''/^e-"^°^^(^-'"). (7) 

dy Stv J J 




Fig. 1(a) Rapidity distribution from a single thermal fire-ball with transverse flow. 
(6) Superposition of uniformly distributed kaon-fire-balls 




Fig. 2 Shematic plot of the fire-balls distributed uniformly 
in longitudinal phase space 
In Eqn.(7) the collective expansions in both the longitudinal and the transverse 
directions have been taken into account. For the convenience of studying the non- 
uniform longitudinal flow, we change the order of integration in Eqn.(7) 

I dr]F{y-r]i), (8) 
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(9) 



The function F{y) plotted in Fig.l can be interpreted as the rapidity distribution from 
a ^^ single thermal fire-hall with transverse flow" and the rapidity distribution dN/dy of 
Eqn's. (7), (8) is the sum of the contributions of a series of such fire-balls with centers 
distributed uniformly in the rapidity range [— r/o, ^70], cf. the shematic plot in Fig. 2 and 
Fig.l(6). 

It is well known that Eqn.(7), or equivalently Eqn's. (8) (9), does not reproduce the 
central dip in the rapidity distirbutions of heavier particles in central relativistic heavy 
ion collisions observed experimentaly already at the AGS energies and confirmed further 
in the experiments at SPS energies. In order to account for this experimental finding 
Ref.p introduced a rapidity dependence of chemical potential /x = ^{rj). However, the 
physical meaning of this dependence is not so clear. 

In the following we will follow the reasoning proposed in Ref.0, i.e. we assume 
that the fire-balls produced in relativistic heavy ion collisions will keep some memory 



on the motion of the incident nuclei, and therefore the distribution of fire-balls, instead 
of being uniform in the longitudinal direction, is more concentrated in the direction 
of motion of the incident nuclei. This means that the distribution of fire-ball is more 
dense at large absolute value of rapidity, as sketched schematically in Fig. 3. 




Fig. 3 Schematic plot of the fire-ball distribution 
in non-uniform longitudinal flow 
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Fig. 4 The dependence of longitudinal distribution of fire-ball 
center p(yie) on ellipticity parameter e 

Using a simple geometrical picture to parametrize the non-uniform longitudinal 
flow 1^, a non-flat distribution function of flre-ball center yie 



piVh 



1 + sinh^(|/ic) 



is introduced and the rapidity distribution (8) becomes 



(10) 



S + S 200 A GeV 




Fig. 5 Rapidity distributions of kaon (a), net proton (b) and lambda (c) 
for central S-S collisions at 200 A GeV/c. Data are taken from Ref.| 
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In Eqn.(lO) e is the "ellipticity" parameter describing the degree of non-uniformity 
of longitudinal flow 0. The dependence oi p{yie) on e is shown in Fig. 4. It can be seen 
from the figure that the larger is the parameter e, the flatter is the distribution p{yic) 
and the more uniform is the longitudinal-flow distribution, cf. Fig. 3. When e —>■ 1, 
the longitudinal-flow distribution is completely uniform (p(yie) — > 1), returning back 
to the cylindrically symmetrical collective flow model, Eqn.(7). 

III. Comparison with experiments 



We present in Fig. 5 (a, b and c) the rapidity distributions of kaon, net proton and 
lambda for central S-S collisions [Q at 200 A GeV/c, respectively. In view of the large 
statistical errors of the experimental data, the fitted values of parameter e are given in 
a region and accordingly the rapidity distributions have a strip shape. 




S "'■.('^) Pb+Pb (158 AGeV/c) 
a 2 - -. 
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Fig. 6 The dependence of longitudinal distribution of fire-ball 
center p{yic) on ellipticity parameter e 

The rapidity limit yieo and the ellipticity e used in the calculation are listed in Table 
I and illustrated in Fig.6 (a). The parameter T is chosen to be 0.12 GeV. 

Table I The values of model-parameters 



Parameter 


S-S Collisions 


Pb-Pb Collisions 


k+ 


p-p 


A 


p-p 


e 


0.7-0.9 


0.3-0.45 


0.7-0.9 


0.56 


yico 


1.98 


1.98 


1.98 


1.96 


< A > 


0.3c 


0.3c 


0.3c 


0.5c 



Since kaons and lambdas are produced through the interaction of colliding nuclei, 
they have less memory on the motion of the incident nuclei. Therefore, the values of 
ellipticity e for kaon and lambda are bigger than that for proton, as shown in Table I. 

It can be seen from the figures that our longitudinally non-uniform collective ex- 
pansion model (NUCEM) reproduces the central dip of rapidity distributions for kaon, 
lambda and net proton very well, in agreement with the experimental findings. 

The appearance or disappearance of central dip is insensitive to the rapidity limit 
j/ieo but depends strongly on the magnitude of the ellipticity e and the mass m of 
produced particles. When transverse flow exists there is a shallow dip also for the 
rapidity distribution of light particle (kaon), cf. Fig.5, which is slightly different from 
the case [0 of one-dimensional longitudinal flow. 

The rapidity distribution of net protons for Pb-Pb collisions at 158 A GeV/c 
is shown in Fig. 7, and the longitudinal distributions of fire-ball center are showm in 
Fig.6(6) for various values of ellipticity e. The solid line in Fig. 7 corresponds to the 



results of our model with ellipticity parameters e = 0.56. The fitted x^/DF is given in 
Fig. 8. The value e = 0.56 corresponds to the minimum of x^- The region of e for x^ 
to increase one unit from the minimum is between 0.40 and 0.745. 



PbtPb158AGeV 





0,7 0.8 0.9 



Fig. 7 Rapidity distribution of net protons 



Fig. 8 The fitting x^ 



Comparing Fig's. 5 and 7 we can see that the net-proton (p— p) rapidity distribution 
is narrower for Pb-Pb |^, than for S-S collisions |^. This indicates an increasing baryon 
stopping for Pb-Pb collisions, giving accordingly a smaller rapidity limit yieo and a 
larger average transverse velocity and uniformity for Pb-Pb interaction. The width of 
the rapidity distributions is mainly controlled by the amplitude yieo of the longitudinal 
flow. For smaller colliding system (S-S collision) a single value of yieo can account for 
the wide distribution of heavier particles (net protons and A) and at the same time fit 
the kaon-distribution well. For the larger colliding system (Pb-Pb), the ?/ieo is smaller 
together with a larger average transverse flow velocity. These are shown also in Table 
I. 



IV. Summary and Conclusions 

Let us begin with a discussion of the thermal freeze-out points. Freeze-out marks 
the transition from a strongly coupled system, which evolves from one state of local 
thermal equilibrium to another, to a weakly coupled one of essentially free-streaming 
particles. If this transition happens quickly enough, the thermal momentum distribu- 
tions (superimposed by collective expansion flow) are frozen in, and the temperature 
and collective flow velocity at the transition point can be extracted from the measured 
momentum spectra. In high energy heavy-ion collisions the freeze-out process is trig- 
gered dynamically by the accelerating transverse expansion and the very rapid growth 
of the mean free paths as a results of the fast dilution of the matter [|1^. Idealizing 
the kinetic freeze-out process by a single point in the phase diagram is therefore not 
an entirely unreasonable procedure. 



In high energy heavy-ion colhsions, due to the transparency of nucleus the partici- 
pants will not lose the historical memory totally and the produced hadrons will carry 
some of their parent's memory of motion, leading to the unequivalence in longitudinal 
and transverse directions, i.e. the flow of produced particles is privileged in the longi- 
tudinal direction. As the increasing of duration time and the pressure of the system 
at SPS energy, it is reasonable to assume that the system expands collectively both in 
the transverse and in the longitudinal directions. This picture has been used by lots 
of models ||12[;1I3]- Fo^' asymptotically high energies, a boost-invariant longitudinal 



expansion model is postulated by Bjorken |]T3[, which gives a plateau in the rapidity 
distribution of produced particles. For finite collision energies such as CERN SPS en- 
ergy or below, a cylindrically symmetrical expansion model was firstly postulated by 
Schnedermann, Sollfrank and Heinz [0 by introducing a cut in rapidity. In this model 
a set of fire-balls with centers located uniformly in the rapidity region [-yieo, 2/ieo]) as 
sketched schematically in Fig. 2, represents the longitudinal flow, and at the same time 
there is radial flow developed simultaneously. It can account for the wider rapidity 
distribution when compared with the prediction of the pure thermal isotropic model 
but failed to reproduce the central dip in the proton and A rapidity distributions. 

In this paper, we argue that the transparency /stopping of relativistic heavy ion colli- 
sions should be taken into account more carefully. It will not only lead to the anisotropy 
in longitudinal-transverse directions, but also render the fire-balls to concentrate more 
in the direction of motion of the incident nuclei. A non-uniform longitudinal flow model 
is proposed, which assumes that the centers of fire-balls are distributed non-uniformly 
in the longitudinal phase space. A parameter e is introduced through a geometrical 
parameterization which can express the non-uniformity of flow in the longitudinal di- 
rection, i.e. the centers of fire-balls of produced particles prefer to accumulate in the 
two extreme rapidity regions (yie ~ ±2/ieo) in the c.m.s. frame of relativistic heavy- 
ion collisions, and accordingly the distribution is diluted in the central rapidity region 
[rje ~ 0). Apart from the non-uniform longitudinal flow, a radial flow in the transverse 
direction is also considered in the present paper. 

It is found that the depth of the central dip depends on the magnitude of the 
parameter e and the mass of produced particles, i.e. the non-uniformity of longitudinal 
flow which is described by the parameter e determines the depth of the central dip for 
produced particles. Comparing with one-dimensional non-uniform longitudinal flow 
model [0, the rapidity distribution of lighter strange particle kaon also shows a dip 
due to the effect of transverse flow. 

Smaller colliding systems distinguish themselves from larger ones not primarily by 
the achieved maximal energy density, but by the occupied collision volume in space 
and time. Compared to S-S collisions, Pb-Pb collisions live longer until thermal freeze- 



out, expand more in the transverse direction, develop more transverse collective flow. 
It is found that transverse collective expansion is expected to be more prominent in 
the heavy Pb-Pb collision system, owing to smaller surface-to-volume ratio, the larger 
flreball volume and the longer duration time of expansion. On the other hand, Pb-Pb 
collisions develop less longitudinal flow, which suggests, together with the larger e, a 
larger stopping in the bigger colliding system. 
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